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Autonomic  materials  are  those  that  automatically  respond  to  a  change  in  environmental  conditions,  such  as 
temperature  or  chemical  composition.  While  such  materials  hold  incredible  potential  for  a  wide  range  of 
uses,  their  implementation  is  limited  by  the  small  number  of  fully-developed  material  systems.  To 
broaden  the  number  of  available  systems,  we  have  developed  a  post-functionalization  technique  where 
a  reactive  Ru  catalyst  ink  is  printed  onto  a  non-responsive  polymer  substrate.  Using  a  succinimide- 
amine  coupling  reaction,  patterns  are  printed  onto  co-polymer  or  biomacromolecular  films  containing 
primary  amine  functionality,  such  as  polyacrylamide  (PAAm)  or  poly-/V-isopropyl  acrylamide  (PNiPAAm) 
copolymerized  with  poly-/V-(3-Aminopropyl)methacrylamide  (PAPMAAm).  When  the  films  are  placed  in 
the  Belousov-Zhabotinsky  (BZ)  solution  medium,  the  reaction  takes  place  only  inside  the  printed  nodes, 
in  comparison  to  alternative  BZ  systems,  where  Ru-containing  monomers  are  copolymerized  with  base 
monomers,  reactive  printing  provides  facile  tuning  of  a  range  of  hydrogel  compositions,  as  well  as 
enabling  the  formation  of  mechanically  robust  composite  monoliths.  The  autonomic  response  of  the 
printed  nodes  is  similar  for  all  matrices  in  the  BZ  solution  concentrations  examined,  where  the  period  of 
oscillation  decreases  in  response  to  increasing  sodium  bromate  or  nitric  acid  concentration.  A 
temperature  increase  reduces  the  period  of  oscillations  and  temperature  gradients  are  shown  to 
function  as  pace-makers,  dictating  the  direction  of  the  autonomic  response  (chemical  waves). 


1  Introduction 

Autonomic  materials  are  those  that  automatically  respond  to  a 
change  in  environmental  conditions.  These  materials  “sense” 
the  environment  as  an  input  signal  and  respond  accordingly 
without  an  external  guide  making  decisions  on  their  behalf. 
Living  organisms  represent  the  most  complex  examples  of 
autonomic  systems.  However,  simpler  autonomous  materials 
also  exist,  where  one  or  more  type  of  response  is  coupled  with 
one  or  more  environmental  condition(s).‘  Depending  on  the 
properties  of  the  autonomous  material,  potential  for  a  wide 
range  of  devices  has  been  discussed,  such  as  energy  storage, 
actuation,  computation,  self-repairing  coatings,  and  encryp¬ 
tion.^’^  Yet,  many  issues  must  be  understood  before  one  can 
design  and  engineer  a  material  with  a  desired  autonomic 
response.  For  example,  what  are  the  relationships  between 
transport,  autonomic  response  and  mechanical  properties  of 
the  material  that  limit  the  size  or  shape  change  caused  by  the 
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environment-response  relationship?  What  is  the  relationship 
between  shape,  pattern  and  time  scale  associated  with  the 
response?  How  can  multiple  responses  be  constructively 
coupled  within  a  single  material?  These  challenges  along  with 
practical  issues  such  as  material  durability  and  handleability 
necessitate  the  further  development  of  new  material  systems 
with  facile  processing  and  composition  control,  as  well  as  the 
continual  development  of  theory  and  modeling  tools. 

The  development  of  a  sub-class  of  autonomic  hydrogels  that 
convert  chemical  energy  into  an  optical  or  mechanical  response 
has  received  substantial  attention  recently.'’^  ^  Such  biocom¬ 
patible  hydrogels  have  the  potential  for  use  in  biomedical 
applications  including  biomimetic  valves  {e.g.  artificial  heart), 
drug  release,  separation  systems,  and  artificial  skin.  Other 
device  applications  include  microfluidics  (actuators)  and 
micromachines  (mass  transport).*  These  polymer-based  systems 
typically  respond  to  variations  in  one  or  more  environmental 
conditions,  such  as  temperature,  pH,  light,  concentration 
gradient,  or  electrical  pulse;  with  a  response  including  changes 
in  size,  shape  or  color.*  The  most  studied  autonomic  hydrogels 
are  based  on  poly(Af-isopropylacrylamide)  (PNIPAAm)  and  the 
Belousov-Zhabotinsky  (BZ)  reaction.*  The  oscillating  BZ  reac¬ 
tion  occurs  in  a  solution  of  sodium  bromate  (SB),  malonic  acid 
(MA),  and  nitric  acid  (NA),  and  is  characterized  by  the  periodic 
oxidation  and  reduction  of  a  metallic  catalyst  (typically  Ru 
between  the  (II)  and  (III)  oxidation  states).”'**  When  the  metal 
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catalyst  is  incorporated  into  PNIPAAm-based  gels,  the  periodic 
oscillation  of  its  oxidation  state  alters  the  hydrophohicity  of  the 
hydrogel.  The  Ru(iii)  state  is  more  hydrophilic  than  the  Ru(ii) 
state  and  therefore  the  PNIPAAm  swells  when  the  catalyst  is 
oxidized. A  corresponding  change  in  size  or  shape  of  the 
hydrogel  occurs  in  phase  with  the  BZ  reaction.*  For  PAAm-hased 
gels,  an  alternate  mechanism  has  been  proposed  where  inter¬ 
chain  interactions  cause  the  gel  to  swell  for  Ru(ii),  and  deswell 
for  Ru(iii).*'‘  Overall,  the  extent  of  swelling  in  these  autonomic 
hydrogels  depends  on  factors  including  the  total  polymer 
concentration,  the  catalyst  concentration,  the  crosslinking 
density,  and  the  temperature  of  the  reaction.* 

Composite  BZ  hydrogels  are  the  next  step  toward  devices. 
These  hydrogels  contain  both  oscillatory  (active)  and  non- 
oscillatoiy  (non-active)  regions  in  the  same  material.^’*^'** 
Communication  between  the  active  nodes  occurs  via  chemical 
and  mechanical  waves  traveling  through  the  inactive  regions. 
Subsequent  oscillatory  symching  of  the  active  regions  leads  to 
emergent  behavior  and  greatly  increases  the  complexity  and 
functionality  of  the  autonomic  response.  Elucidating  the  design 
rules  associated  with  building  these  composite  hydrogels  has 
been  the  focus  of  recent  experiments  and  simulations.*’*^’*®  “ 
The  details  of  internode  communication  depends  on  the 
properties  of  the  hydrogel  and  the  geometry  of  node  arrange¬ 
ment.*®  Photo  illumination  can  also  be  used  to  modulate 
internode  activity  and  overall  composite  behavior.**^^*  For 
example,  inter-band  transitions  within  Ru  from  absorption  of 
436  nm  light  inhibit  BZ  oscillations  and  enables  continuous 
grey  scale  tuning  of  the  reaction  and  the  hydrogel  response.  By 
combining  these  concepts,  wormlike  motion  and  rudimentary 
swarming  have  been  predicted.^*’^^ 

Unfortunately,  the  materials  available  to  construct 
composite  hydrogels  with  requisite  mechanical  robustness, 
simple  fabrication,  and  precision  of  node  geometry  are 
limited. *’**’“  The  majority  of  BZ  gel  demonstrations  in  the 
literature  are  based  on  PNIPAAm  as  a  monolithic  cube,  film, 
sphere,  or  cylinder  of  tens  of  microns  to  a  few  millimeters  in 
size.*  Recent  approaches  to  overcome  these  limitations  include 
using  base  pol3miers  with  more  robust  mechanical  properties 
such  as  polyacrylamide  (PAAm)**’**  and  gelatin,*’*®  and 
employing  novel  synthesis  techniques  and  catalyst  binding 
chemistries  such  as  UV  mask  patterning,**’**  and  post¬ 
functionalization.*’*®  Of  special  note  with  regard  to  complex 
fabrication  is  the  development  of  autonomic  hydrogels 
compatible  with  additive  processing  (manufacturing)  tech¬ 
niques.  For  example,  the  thermal  gelation  characteristics  of 
gelatin  enable  printing  of  self-supporting  structures  by 
extruding  hot  gelatin  (low  viscosity)  onto  a  cold  substrate  (high 
viscosity).*’*®  Glutaraldehyde  crosslinking  of  co-printed,  or 
printed  and  backfilled,  structures  of  active  (Ru-containing)  and 
non-active  gelatin  result  in  monolithic  systems.  Printing  ther¬ 
mogels  can  be  challenging  however  due  to  their  temperature 
sensitive  viscosity.  Current  reports  are  limited  to  ~100  pm 
resolution,  which  may  be  inadequate  for  certain  applications.*® 
To  further  expand  the  processing-material  concepts  for 
autonomic  hydrogel  fabrication,  we  demonstrate  herein  reac¬ 
tive  printing  of  Ru-containing  BZ-hydrogel  composites.  Using  a 
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succinimide-amine  coupling  reaction,  responsive  Ru  patterns 
are  embedded  within  co-polymer  films  containing  primary 
amine  functionality,  such  as  PAAm  or  PNIPAAm  copolymerized 
with  Af-(3-aminopropyl)methaciylamide  (APMAAm).  Reactive 
printing  reduces  the  number  of  steps  to  create  the  active- 
inactive  node  patterns  to  two;  one  polymerization  step  to  create 
a  cross-linked  monolithic  film,  and  a  second  printing  step  to 
create  the  reactive  patterns.  The  Ru  catalyst  is  delivered  in  an 
ink  solution  of  DMF,  which  can  be  applied  with  something  as 
simple  as  a  micro-dispenser  or  inkjet  printer. 

2  Experimental 

2.1  Materials 

iV-Isopropyl  acrylamide  (NIPAAm),  acrylamide  (AAm),  methyl- 
enebisacrylamide  (MBA),  ammonium  persulfate  (APS),  tetrame- 
thylethylenediamine  (TEMED),  Ar,JV-dimethylformamide  (DMF), 
and  bis(2,2'-bipyridine)-4'-methyl-4-carboxybipyridine-ruthe- 
nium  N-succinimidyl  ester-bis(hexafluorophosphate)  (Ru(sbpy)) 
were  purchased  from  Sigma-Aldrich  and  used  as  received.  Af-(3- 
Aminopropyl)methaciylamide  hydrochloride  (AMPAAm)  was 
purchased  from  Polyscience,  Inc.  and  used  as  received.  For  BZ 
reactants,  sodium  bromate  (SB)  and  malonic  acid  (MA)  were 
purchased  from  Sigma-Aldrich  and  used  as  received.  Nitric  acid 
(NA)  was  purchased  from  Sigma-Aldrich  as  70%  reagent  and 
diluted  to  the  concentration  used  in  the  reactions. 

2.2  Reactive  printing 

Substrate  polymers  are  fabricated  into  thin  films  ~400  pm  thick 
by  the  following  example  procedure.  100  mg  of  AAm,  60  mg 
APMAAm,  3.2  mg  MBA  crosslinker,  5  mg  APS,  and  10  pL  of 
TEMED  were  dissolved  in  1  mL  of  deionized  (DI)  water.  Other 
variations  of  this  procedure  include  using  NIPAAm  in  place  of 
AAm,  using  a  90  :  10  mixture  of  NIPAAm-AAm,  and  adjusting 
the  mass  of  MBA  cross-linker,  for  example,  to  1.6  mg.  The 
solution  was  injected  between  two  glass  panels  covered  with 
hydrophobic  BYTAC  and  separated  by  4  BYTAC  layers 
(~400  pm).  This  thickness  provided  adequate  mechanical 
rigidity  for  subsequent  handling  and  processing.  The  glass 
panels  are  clamped  together  and  the  film  is  allowed  to  poly¬ 
merize  for  4  hours  at  room  temperature.  After  pol3nnerizing,  the 
film  is  allowed  to  dry  for  20  minutes  and  cut  with  a  razor  blade 
into  square  sections  approximately  2  cm  x  2  cm.  These  square 
sections  contain  approximately  90%  water  by  weight  when  fully 
hydrated.  1.5  mg  of  Ru(sbpy)  catalyst  is  dissolved  in  30  pL  of 
DMF  to  make  the  reactive  ink  solution.  For  experiments  to 
determine  the  period  dependence  on  temperature  and  reactant 
concentration,  one  drop  (approx.  5  pL)  of  reactive  ink  is  applied 
by  a  microliter  scale  fluid  dispenser  equipped  with  a  26.5  gauge 
needle  to  the  square  section  of  the  substrate  to  form  a  circular 
spot  in  the  center  with  diameter  ~l-2  mm.  To  demonstrate  the 
ability  to  create  complex  shapes  and  patterns,  we  also  printed 
spots  of  various  sizes,  lines,  letters,  and  geometrical  shapes 
such  as  squares,  triangles  and  circles.  The  procedure  used  to 
print  these  shapes  is  the  same  as  the  1-2  mm  spots  above,  but 
the  number  of  drops  is  either  increased  (for  large  spots)  or  the 
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drops  overlapped  each  other  to  form  a  line  (for  lines  and 
shapes).  For  line  printing,  a  flow  rate  of  5  pL  s^^  and  a  dispenser 
speed  of  2  mm  s^^  is  used.  After  application  of  the  ink,  the 
substrates  are  placed  in  an  oven  at  35  °C  and  covered  for  3  hours 
to  allow  the  succinimide-amine  coupling  reaction  to  complete. 
By  3  hours,  the  films  have  dried  to  less  than  5%  water  and  the 
coupling  reaction  has  completed.  After  this,  they  are  rehydrated 
and  washed  with  DI  water  to  remove  unreacted  monomer  and 
excess  Ru(sbpy)  until  the  wash  no  longer  turns  orange  and 
remains  clear.  The  washing  process  takes  5-10  cycles  of 
immersing  in  ~10  mL  of  fresh  water  for  10  minutes  followed  by 
removal  of  the  contaminated  rinse. 

2.3  Belousov-Zhabotinslg?  (BZ)  driven  oscillations 

The  printed  sections  are  placed  in  a  small  Petri  dish  with  5  mL 
of  BZ  reactant  solution  and  covered  to  prevent  solvent  evapo¬ 
ration.  We  note  that  failure  to  cover  the  reaction  vessel  results  in 
solvent  evaporation,  and  will  lead  to  increasing  reactant 
concentrations  over  time  and  an  apparent  long-time  transitory 
behavior  (time-dependent  period  of  oscillations  and  shrinking 
of  the  hydrogel).  In  the  BZ  reaction,  an  organic  compound  (MA) 
is  oxidized  by  an  inorganic  oxidant  (SB),  in  the  presence  of  a 
strong  acid  (NA)  and  a  catalyst  (Ru).  BZ  reactant  concentrations 
range  from  0.08-0.25  M  for  SB,  0.04-0.14  M  for  MA,  and  0. 5-1.1 
M  for  NA  with  the  base  concentrations  being  0.2  M  (SB),  0.1  M 
(MA)  and  1.0  M  (NA).  One  reactant  is  varied  at  a  time  while  the 
others  are  kept  at  the  base  concentration  in  order  to  determine 
the  effect  of  each  individual  reactant  on  the  period  of  chemical 
waves.  Once  in  the  BZ  solution,  the  color  of  the  spots  changes 
from  orange  to  transparent  as  the  Ru  oscillates  between  the  (II) 
and  (III)  oxidation  states,  respectively.  The  BZ  reaction  is 
observed  as  swirls  or  lines  of  color  that  transverse  the  spots  at 
regular  intervals  and  are  captured  by  images  taken  with  a  Dino- 
Lite  microscope  camera  spaced  2  s  apart.  Total  image  acquisi¬ 
tion  and  reaction  time  is  5  hours.  The  temperature  of  the 
reaction  is  controlled  with  a  Peltier  heating  plate  under  the  Petri 
dish  connected  to  a  variable  voltage  DC  power  source  and 
monitored  by  a  thermocouple  adhered  to  the  plate.  The  Dino- 
Lite  camera  images  are  separated  into  red,  green  and  blue 
channel  8-bit  images  and  the  intensity  of  a  single  pixel  of  the 
blue  channel  is  analyzed  in  Matlab  to  find  the  frequency  of 
the  color  oscillations.  The  exact  location  of  the  pixel  chosen  on 
the  printed  spot  had  no  effect  on  the  results. 

3  Results  and  discussion 

3.1  Reactive  printing 

A  diagram  demonstrating  the  ink  printing  procedure  and  the 
succinimide-amine  coupling  reaction  between  Ru(sbpy)  and 
APMAAm  units  is  shown  in  Fig.  1.  In  brief,  the  succinimide- 
amine  reaction  forms  an  amide  bond  between  the  nitrogen 
atom  in  the  amine  group  on  the  co-polymer  and  the  acyl  residue 
on  the  Ru  complex.  The  reaction  takes  place  in  water  and  is 
most  efficient  at  pH  ~7-8  and  temperatures  ~40  A  side 

hydrolysis  reaction  of  the  ester  group  competes  with  the  amide 
bond  formation  in  water  which  lowers  the  conjugation  yield. 
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Fig.  1  (Above)  Schematic  diagram  of  the  reactive  printing  procedure. 
Patterns  can  be  programed  by  an  inkjet  printer,  and  complex  patterns 
such  as  spots,  lines,  shape  outlines,  and  letters  have  been  printed. 
(Below)  Diagram  showing  the  succinimide-amine  coupling  reaction 
adapted  from  ref.  27.  The  succinimide-ester  functional  group  on  the 
Ru  complex  reacts  with  the  amine  groups  of  the  APMAAm  units  along 
the  co-polymer,  which  forms  an  amide  bond  that  couples  the  Ru  to 
the  matrix.  The  structure  of  APMAAm  is  also  shown. 

Herein  a  microliter  dispenser  with  a  26.5  gauge  needle  was  used 
to  meter  5  pL  drops  to  produce  plxelated  patterns.  Uniform 
lines  were  created  by  correlating  the  rate  of  fluid  dispensing 
(5  pL  s^’^)  with  lateral  translation  (2  mm  s^^).  The  process  could 
be  easily  Integrated  Into  a  programmable  Inkjet  printer  allow¬ 
ing  a  nearly  limitless  number  of  shapes  and  sizes  of  reactive 
patterns.^®  Examples  of  shapes  and  letters  printed  with  this 
procedure  are  shown  in  Fig.  2,  and  a  movie  showing  how  the 
printed  “letters”  behave  in  the  BZ  reaction  is  found  in  the  ESI, 
Movie  Sl.f  Four  key  factors  determine  the  response  of  these 
printed  composite  hydrogel  patterns  In  the  BZ  solution.  As 


Fig.  2  Examples  of  shapes  printed  onto  PAAm-co-APMAAm 
substrates  using  the  procedures  described  in  the  text.  Dashed  black 
outlines  have  been  drawn  on  the  images  to  show  the  edge  of  the 
substrate  material  in  solution  and  to  clarify  film  borders  due  to 
reflections.  Solid  orange  lines  have  been  drawn  on  the  circle  and 
square  printed  shapes  to  help  show  where  the  ink  was  placed. 
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discussed  in  the  following  sections,  these  include  the  depth  or 
penetration  of  the  printed  features,  the  concentrations  of  BZ 
reactants,  the  mechanical  properties  of  the  substrate  material, 
and  the  temperature  of  the  reaction. 

3.2  Feature  depth 

The  lateral  size,  and  most  importantly  depth,  of  the  printed 
active  node  depends  on  the  extent  to  which  the  catalyst-con¬ 
taining  ink  diffuses  into  the  substrate  material.  The  profile  will 
influence  reaction  kinetics  since  reactant  diffusion  rates  are 
slower  inside  the  polymer  matrix  than  at  the  surface.  It  also 
determines  whether  the  chemical  waves  are  occurring  at  the 
film  surface  or  in  the  bulk,  and  strongly  influences 
the  mechanical  response  of  the  active-inactive  patterns,  as  the 
inactive  volume  below  the  active  node  will  act  as  a  constraining 
boundary  to  the  autonomic  swell-deswell  behavior.  In  general, 
the  profile  of  the  printed  node  is  directly  related  to  the  mutual 
solubility  of  the  reactive  Ru(sbpy),  the  ink,  and  the  slightly 
swollen  hydrogel  substrate.  If  the  polymer  and  Ru(sbpy)  are 
both  soluble  in  the  solvent  used  for  the  ink,  the  printed  feature 
will  quickly  spread,  but  also  penetrate  deeply  into  the  hydrogel 
network.  This  solubility  combination  will  maximize  the  total 
Ru(sbpy)  incorporated  into  the  hydrogel  before  hydrolysis 
deactivates  the  coupling  reaction  between  pendent  amine  and 
Ru(sbpy),  but  will  compromise  the  ability  to  print  well-defined 
features.  As  an  optimization  trade,  the  inks  formulated  for  Fig.  2 
used  DMF,  a  solvent  in  which  PAAm  and  PNIPAAm  has  low 
solubility.^*  The  edges  of  these  printed  features  are  crisp  and 
well-defined.  Flowever,  since  the  DMF  solution  does  not  readily 
diffuse  into  the  polymer,  the  coupling  reaction  only  takes  place 
near  the  surface.  The  reactive  ink  penetrates  only  to  a  depth  of 
10-20  |im,  as  demonstrated  in  the  cross  section  of  a  printed 
active  node  within  PAAm-co-APMAAm  shown  in  Fig.  3. 
Depending  on  the  intended  application,  it  is  likely  that  the 
reactive  ink  can  be  formulated  to  be  optimized  for  a  given 
substrate  and  the  required  lateral  and  depth  resolution  of  the 
active  node  by  modifying  co-solubility  {e.g.  DMF  :  water  ratio  of 
ink),^*  Ru(sbpy)  concentration  in  the  ink,^’“  polymer 


Fig.  3  Profile  of  reactively  printed  active  node.  An  active  node  in 
PAAm-co-APMAAm  substrate  material  was  cut  into  a  strip  and  optical 
microscope  images  were  taken  of  the  top  and  side  faces.  The  cross 
sectional  side  face  image  shows  the  penetration  depth  of  the  ink  of 
approximately  10-20  |im  and  the  top  face  shows  the  well-defined 
border  of  the  printed  node. 


concentration  in  the  substrate  (water  content  of  film),  and  the 
ratio  of  APMAAm  to  base  polymer  in  the  hydrogel. 

3.3  Oscillation  behavior  of  composite  BZ  hydrogel 

The  oscillation  behavior  of  the  printed  composite  hydrogel 
depends  on  the  reactant  concentration  of  the  BZ  solution,  the 
mechanical  characteristics  of  the  composite,  and  temperature. ‘ 
The  dependence  on  reactant  concentration  and  temperature 
provide  potential  means  to  control  the  period  of  oscillation 
across  the  nodal  array  by  making  uniform  changes  to  these 
inputs.  Using  gradients,  such  as  with  temperature,  provides 
control  of  the  direction  of  the  chemical  waves. 

Fig.  4  summarizes  the  dependence  of  the  oscillatory 
behavior  at  20  °C  of  single  active  nodes  on  the  concentration  of 
BZ  chemical  reactants.  The  reactant  concentrations  were  varied 
individually  while  holding  the  other  two  reactants  constant  at 
the  base  concentration  (see  Experimental  section).  The  rela¬ 
tionship  between  the  period  and  reactant  concentration  is  given 
by  the  equation  P  =  A[SB]*[MA]'^[NA]‘^  where  P  is  the  period,  A  is 
a  reaction  coefficient  and  b,  c,  and  d  are  exponents  independent 
of  the  other  reactant  concentrations.  As  shown  by  Smith  and 
coworkers  with  gelatin^  and  Yoshida  and  coworkers  with  PNI- 
PAAm,‘*’*“  the  response  empirically  collapses  to  a  single  plot  by 
using  a  sum  of  logs,  P  =  b  ln[SB]  +  c  ln[MA]  +  d  ln[NA].  Overall, 
the  period  shortens  with  increasing  reactant  concentrations. 
However,  the  period  has  a  much  stronger  dependence  on  SB 
and  NA  than  MA.  The  series  of  images  at  the  top  of  Fig.  4  show  a 


Fig.  4  Period  of  oscillations  at  20  “C  of  printed  active  nodes  (2-5  mm) 
in  various  hydrogels.  Each  hydrogel  contains  38%  by  weight  of 
APMAAm  co-monomer.  The  concentrations  of  the  BZ  reactants  are 
varied  from  0.08-0.25  M  for  [SB],  0.04-0.14  M  for  [MA],  and  0.5-1.1  M 
for  [NA].  The  solid  line  represents  the  nonlinear  least  squares  fit  to  the 
period  based  on  the  reactant  concentrations  expressed  as  a  sum  of 
logs.  The  dashed  line  represents  the  behavior  of  the  BZ  reaction  in  the 
absence  of  hydrogel  and  with  the  same  concentration  of  unbound  Ru 
catalyst.’*  The  series  of  images  in  the  inset  show  an  example  of  a 
chemical  wave  traveling  across  a  printed  node.  Error  bars  indicate 
variation  of  the  period  within  one  experiment. 
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representative  single  wave  response  moving  across  an  active 
node  of  PAAm-co-APMAAm.  The  color  has  heen  converted  to 
greyscale  to  enhance  the  contrast  between  the  (II)  state,  In 
hlack,  and  (III)  state,  in  white.  An  example  of  the  color  oscilla¬ 
tions  over  time  of  a  printed  spot  are  shown  in  the  ESI,  Fig.  Slf 
and  a  movie  showing  an  example  of  the  behavior  for  a  spot  is 
shown  in  Movie  S2.t 

No  significant  difference  in  oscillatory  behavior  was  found 
between  reactively  printed  nodes  in  PAAm-co-APMAAm,  PNI- 
PAAm-co-APMAAm,  and  a  10  :  90  copolymer  of  PAAm  and  PNI- 
PAAm.  Each  hydrogel  contains  38%  by  weight  of  APMAAm. 
Repeated  experiments  indicate  that  sources  of  non-negligible 
experimental  scatter  include  fluctuations  in  lab  temperature, 
lab  humidity,  and  reactant  concentrations.  These  occur  during 
film  synthesis,  node  printing,  and  preparation  and  storage  of 
the  BZ  solution.  However,  the  largest  source  of  noise,  as  indi¬ 
cated  by  the  error  bars  in  Fig.  4,  is  the  drift  of  the  period  length 
from  cycle  to  cycle  within  a  single  experiment.  Considering 
these  points,  the  results  with  different  substrate  materials  is 
remarkably  consistent.  Additionally,  compared  to  the  reported 
period  of  oscillation  of  free  Ru  catalyst  within  a  BZ  solution  at 
the  same  temperature  (dashed  line.  Fig.  4),  there  is  again  little 
difference  with  the  overall  trend  of  the  data.“’“  This  implies 
that  although  the  catalyst  is  chemically  bound  to  the  matrix,  the 
reaction  kinetics  are  no  different  than  those  found  in  a  solution 
with  unbound  catalyst.  Similar  results  for  the  period  of  oscil¬ 
lations  were  determined  for  lateral  node  diameters  from  1  mm 
to  5  mm  as  well  as  for  shapes  and  lines,  even  though  the  shape 
of  the  wave  differed  with  the  size  and  shape  of  the  node  {i.e.  one 
directional  waves  versus  spiral  waves). 

Elucidating  the  impact  of  matrix  composition  on  oscillatoiy 
response  is  not  straightforward,  as  crosslink  and  hydrogel 
concentration  has  also  been  shown  to  play  a  role.  For  example, 
previous  studies  with  gelatin^  demonstrated  a  different  expo¬ 
nential  dependence  of  the  period  on  the  reactant  concentra¬ 
tions.  The  exponents  for  SB,  MA  and  NA  in  gelatin  were  found  to 
be  —0.512,  1.387,  and  —0.348  for  a  400  pm  cube,  which 
compares  to  1.085,  0.278,  and  1.074  for  this  work.  At  the  pH  of 
the  BZ  solution  (1-2),  gelatin  is  a  charged  polyelectroly^te  due  to 
ionization  of  amino  acids  residuals  such  as  glycine,  glutamic 
acid  and  alanine.  Also,  physical  gelation  of  the  gelatin  creates  a 
molecularly  inhomogeneous  network.  Both  factors  could  result 
in  different  diffusion  characteristics  for  the  BZ  reactants. 
Extrapolating  these  prior  observations  to  similar  BZ  reactant 
concentrations  examined  herein  would  predict  larger  periods  in 
the  range  of  60-300  s.  Similarly,  reports  by  Yoshlda  and 
coworkers,  using  monolithic  PNIPAAm  gels  with  similar  cross- 
linking  and  polymer  concentrations,  would  imply  long  periods 
of  100-500  However,  experiments  by  Nuzzo  and  coworkers 
with  PAAm  monoliths  with  much  higher  crosslinking  (~12%  by 
mass)  than  discussed  here,  reported  similar  period  lengths  in 
the  range  of  10-60  s.*"* 

As  such,  there  are  two  potential  explanations  for  the  simi¬ 
larity  in  chemical  wave  period  shown  in  Fig.  4.  One  is  that  for 
the  printed  composite  hydrogels,  the  reaction  takes  place  only 
in  the  near  surface  region  due  to  the  limited  penetration  of  the 
active  node.  Because  of  this,  possible  reduction  of  diffusion 
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rates  of  the  BZ  reactants  by  the  surrounding  inactive  hydrogel  is 
not  a  factor.  The  comparable  oscillation  rate  with  free  solution 
bolsters  this  explanation.  Thus  the  “matrix  independent 
response"  is  a  geometrical  effect.  Alternatively,  the  relatively  low 
matrix  concentration  (10%  polymer  by  volume  when  hydrated) 
and  neutrality  of  PAAm  does  not  substantially  hinder  the 
diffusion  of  the  BZ  reactants  relative  to  the  surrounding 
aqueous  solution.  Prior  permeability  studies  of  interpenetrating 
PAAm-PNIPAAm  hydrogels  show  that  rates  of  water  uptake  and 
diffusion  of  tracer  Orange  II  azo  dye  depend  non-monotonically 
on  temperature,  crosslinking,  and  total  polymer  content.^^ 
However,  at  temperatures  below  the  LOST  (32  °C),  gel  perme¬ 
ability  of  swollen  films  of  comparable  polymer  and  crosslinking 
concentration  used  herein  only  displayed  a  small  effect  on  the 
relative  PNIPAAm  composition.  This  implied  that  the  impact  of 
the  matrix  composition  was  secondary  to  the  absolute  polymer 
concentration  in  these  formulations.  Further  studies  using 
print-based  fabrication  techniques  are  required  to  sort  out  the 
interplay  between  the  relative  size  and  shape  of  the  active  node 
on  the  oscillatory  behavior  of  the  composite. 

The  reactively  printed  composite  hydrogels  displayed  oscil¬ 
latory  strains  up  to  0.5%.  This  was  irrespective  of  the  extent  of 
PNIPAAm  incorporation  into  the  PAAm  polymer  backbone.  This 
is  in  contrast  to  reports  of  monolithic  gelatin,^  PNIPAAm,‘’*"’“ 
and  PAAm“  based  autonomic  hydrogels,  with  oscillatory 
swelling  range  from  2-20%.  As  noted  above,  the  limited  pene¬ 
trability  of  the  Ru  ink  that  results  in  the  crisp  pattern  formation 
also  hinders  the  ability  to  uniformly  distribute  the  Ru  catalyst 
throughout  a  monolith  of  sufficient  thickness  and  robustness  to 
handle.  Attempts  to  apply  ink  to  all  sides  of  a  cube  (400  pm^) 
resulted  in  a  core-shell  structure  and  no  measurable  swelling. 
As  discussed,  the  printing  depth  of  the  Ru  catalyst  is  inversely 
related  to  the  precision  of  the  print  geometry.  This  is  also 
effectively  a  trade-off  between  print  precision  and  swell  perfor¬ 
mance  of  the  BZ  gel.  For  comparison  purposes  we  normalize 
print  depth,  h,  by  gel  thickness,  H,  to  define  fractional  pene¬ 
tration  depth  hlH.  The  print  depth  and  gel  thickness  in  our 
system  are  h  =  10-20  pm  and  H  =  400  pm  respectively,  resulting 
in  penetration  depths  of  less  than  0.05. 

Fig.  5  summarizes  finite  element  (FE)  simulations  (see  SI 
methods)  to  provide  a  better  understanding  of  the  bulk  swelling 
behavior  of  notional  BZ  gel  cubes  with  variable  penetration 
depths  of  active  isotropic  material.  As  shown  in  Fig.  5A,  we  see 
that  for  a  0.05  penetration  depth  the  swell  ratio  of  the  entire  gel 
Is  low,  even  if  the  active  material  swells  10%.  For  an  elastically 
incompressible  material,  bulk  swelling  performance  is  essen¬ 
tially  a  volume-weighted  average  of  the  active  and  non-active 
components.  As  the  penetration  depth  of  the  active  node 
decreases,  the  influence  of  the  non-active  underlying  region  is 
more  pronounced.  Bulk  swelling  could  be  Increased  by  using  a 
less  stiff  non-active  gel.  This  however,  will  affect  the  overall 
mechanical  robustness  of  the  composite  structure,  potentially 
weakening  it  beyond  practical  transport  and  positioning 
requirements. 

Fig.  5  also  indicates  that  the  swelling  of  the  composite 
hydrogel  Is  heterogeneous  due  to  the  constraint  of  the  non- 
active  material.  This  heterogeneity  is  distinctly  seen  in  the  drop 
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Fig.  5  Mechanical  constraints  on  autonomic  swell-deswell  behavior 
of  reactively  printed  node.  (A).  Bulk  swell  behavior  transitions  from 
constrained  to  free  swelling  as  the  active  node  penetrates  the  gel 
thickness,  H.  Insets  show  swollen  configuration  of  varying  penetration 
geometries.  (B)  Plot  of  in-plane  strain  near  the  surface  of  the  active 
material  versus  penetration  depth.  In-plane  swelling  (e^x,  CzJ  near  the 
surface  of  the  active  material  is  constrained  at  low  penetration  depths. 
The  black  dot  indicates  where  the  strain  values  were  taken.  Figures  to 
the  right  are  contour  plots  of  exx  for  different  h/H  ratios  showing  strain 
distribution  (uniform  color  gradient  between  red:  10%  and  blue:  <0%). 


of  in-plane  strains  in  the  center  of  the  cube  near  the 

active  material  surface  (Fig.  5B).  The  proximity  of  this  top 
surface  to  the  active-non-active  interface  at  low  penetration 
depths  limits  the  in-plane  strain  of  the  active  material.  The 
strain  normal  to  the  surface  (sj^,)  compensates  for  this 
constraint  by  swelling  beyond  the  prescribed  10%,  resulting  in  a 
mushroom  geometry.  As  seen  in  the  contour  plot  of  Fig.  5B,  the 
uniformity  of  the  in-plane  strain  field  in  the  active  material 
increases  as  the  penetration  depth  increases  and  the  active- 
non-active  interface  moves  farther  away  from  the  surface.  Thus 
at  short  penetration  depths,  oscillation  of  surface  “bumps”  in 
the  center  of  the  active  nodes  is  more  likely  than  overall 
macroscopic  strain  of  the  composite  hydrogel. 

Recent  studies  Indicate  that  mechanical  forces  can  resusci¬ 
tate  chemical  oscillations  in  BZ  systems^  and  swelling  may  be 
able  to  coordinate  tbe  chemical  oscillations  of  neighboring  Ru 
nodes. Mechanical  triggering  of  neighboring  nodes  will  likely 
require  a  minimal  in-plane  strain  magnitude,  which  could  be 
masked  at  low  penetration  depths  by  the  non-active  material. 
Compensatory  swelling  normal  to  the  surface  (£j,j,),  will  be 
ineffective  at  stimulating  S3mch  as  its  action  is  orthogonal  to  the 
plane  of  active  nodes.  Thus  for  optimal  mechanical  syncing  of 
planar  arrays,  the  active  nodes  should  span  the  film  thickness. 
Chemical  synching  of  active  nodes,  which  depends  on  reactant 
diffusivity,  will  be  influenced  more  by  relative  placement  of 
nodes  than  penetration  depth.  These  trade-offs  highlight  the 
importance  of  clearly  identifying  the  design  objectives  of  a 
composite  BZ  device  before  selection  of  a  fabrication  technique 
and  material. 

Finally,  control  of  the  oscillation  behavior  of  the  reactively 
printed  nodes,  including  both  period  and  wave  direction,  is 
demonstrated  in  Fig.  6  and  7  through  the  use  of  a  temperature 
jump  or  gradient.  As  an  example.  Fig.  6  summarizes  the 
reversible  oscillatory  control  of  an  active  node  in  PAAm-co- 
APMAAm  using  a  simple  Peltier  heating  plate  under  the  reac¬ 
tion  dish  ([SB]  =  0.2  M,  [MA]  =  0.1  M,  and  [NA]  =  1.0  M). 
Changes  in  the  frequency  of  oscillations  were  observed  witbin  a 
few  minutes  of  changing  the  temperature  and  were  reversible. 
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Fig.  6  Oscillatory  control  of  active  node  in  PAAm-co-APMAAm  via 
temperature  jump.  The  period  of  chemical  waves  reversibly  changes 
with  temperatures.  Approximately  25  s  periods  are  seen  at  25  “C.  At 
15  °C  the  period  lengthens  to  ~55  s,  then  returns  to  20-25  s  once  a 
25  °C  temperature  is  restored.  Finally,  the  period  returns  again  to  ''--55  s 
periods  at  15  °C.  Inset  shows  the  active  node  denoting  the  region 
(black  line)  in  which  the  temporal  behavior  is  displayed  above  the  data. 
Node  diameter  is  2  mm  and  dark  regions  indicate  (II)  state  while  light 
regions  indicate  (III). 


Specifically,  30  minutes  after  beginning  the  reaction,  the 
temperature  of  the  plate  was  changed  from  25  °C  to  15  °C.  Then 
at  140  minutes  from  15  °C  to  25  °C,  and  finally  at  200  minutes 
back  to  15  °C.  Within  approximately  15  minutes  of  the 
temperature  changes,  the  reaction  reaches  a  new  steady  state 
with  an  oscillation  period  of  55  seconds  (15  °C)  or  25  seconds 
(25  °Cj.  For  these  reactant  concentrations  and  for  temperatures 
between  15  °C  and  35  °C,  the  period  of  oscillations  was  found  to 
empirically  follow  an  exponential  dependence  of  the  form 


Fig.  7  Directional  control  of  chemical  wave  in  P/V\m-co-APMAAm  gel 
with  temperature  gradient.  Plot  of  the  relative  phase  lead  and  phase  lag 
of  Ru  oxidation  between  the  black  center  node  and  the  outer  green 
and  blue  nodes,  which  shows  the  transition  from  the  initially  longer 
right-to-left  20  °C  waves  to  the  shorter,  edge-initiated,  left-to-right 
35  °C  waves.  The  left  side  temperature  was  increased  to  35  “C  at 
nominal  time  equals  zero.  Phase  is  normalized  to  period  length, 
P  (P  ~  35  s  at  20  °C).  Schematic  contains  gel  geometry  and  orientation 
of  thermal  gradient.  Grayscale  images  show  signal  intensity  along  the 
yellow  dashed  line  from  (a)  and  (b),  plotted  in  5  minute  sections. 
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P  =  5.7  exp(30(l/r)),  where  P  =  period  in  seconds,  and  T  = 
temperature  in  °C.  Recall  PAAm  does  not  exhibit  a  LCST  such  as 
seen  in  PNIPAAm.®^  Thus  the  temperature  dependence  is 
expected  to  follow  BZ  reaction  behavior,  which  is  typically 
Arrhenius-like.'  Note  that  empirically  similar  temperature 
behavior  has  been  seen  for  actuators  made  from  PNIPAAm; 
however  quick-responding  reversibility  was  not  demonstrated, 
which  is  most  likely  due  to  the  hysteresis  in  swell-deswell 
behavior  of  the  gel." 

The  temperature  dependence  on  the  period  can  also  be  used 
to  control  the  direction  of  the  waves,  as  summarized  in  Fig.  7.  A 
sample  with  a  straight  line  was  placed  in  the  BZ  solution  with 
the  Peltier  plate  set  to  a  uniform  20  °C.  At  this  uniform 
temperature  there  is  no  preference  for  the  direction  of  the  waves 
at  steady  state.  In  Fig.  7,  the  waves  moved  from  right  to  left 
across  the  printed  line.  After  reaching  steady  state,  the  left  half 
of  the  film  was  placed  onto  a  35  °C  plate,  while  the  right  half 
remained  on  the  20  °C  plate.  After  an  initial  equilibration  time, 
waves  began  initiating  at  the  far  left  edge  of  the  sample  and 
moving  from  left  to  right  at  a  shorter  period,  consistent  with  the 
temperature  jump  behavior  in  Fig.  6.  Eventually  the  short  left- 
to-right  waves  took  over  the  longer  right-to-left  waves  along  the 
entire  sample.  The  complete  movie  of  this  behavior  is  included 
in  the  ESI  (Movie  S3t).  This  transition  in  wave  direction  was 
quantified  using  the  relative  phases  between  the  oxidation 
peaks  measured  at  three  uniformly  spaced  locations,  denoted  as 
the  colored  dots  in  the  upper  right  of  Fig.  7.  Initially,  the  wave 
originates  from  the  right  side  of  the  sample  (green  dot).  Then 
from  approximately  25  minutes  to  32  minutes,  BZ  waves 
approach  from  both  sides  of  the  black  node.  Eventually  the  left 
(blue)  side  waves  overcome  the  right  (green)  side  and  lead  the 
wave  direction.  In  other  words,  the  high  temperature  (short 
period)  side  of  the  sample  became  a  pacemaker  for  the  rest  of 
the  sample.  This  is  consistent  with  previous  work  that  showed 
when  nodes  with  two  different  frequencies  are  close  enough  to 
communicate  (and  in  this  case  have  no  separation  between  the 
35  °C  side  and  the  20  °C  side),  the  one  with  the  shorter 
frequency  acts  as  the  pacemaker.'^''^  Further  explanation  of 
wave  direction  analysis  is  given  in  the  ESI,  including  the 
Fig.  S2.t  Extending  this  into  a  2D  or  3D  system  could  provide  a 
unique  means  to  convert  multi-dimensional  thermal  informa¬ 
tion  into  an  optical  signal,  with  potential  application  as  a 
thermocouple  or  anti-tampering  sensor. 

4  Conclusion 

In  this  study,  reactive  printing  is  demonstrated  as  a  facile 
means  to  create  autonomic  composite  hydrogels.  Using  succi- 
nimide-amine  coupling  reaction  and  ink  formation,  fine 
detailed  patterns  of  Ru(sbpy)-containing  active  regions  are 
formed  on  substrates  containing  primary  amines,  such  as 
PAAm  and  PNIPAAm  hydrogels  copolymerized  with  PAPMAAm. 
Importantly,  we  find  that  the  period  of  the  printed  active  node  is 
not  affected  by  the  choice  of  substrate  polymer  and  is  not 
measurably  different  than  that  of  the  Ru  catalyst  in  solution. 
This  demonstrates  that  for  these  composite  systems  and  poly¬ 
mer  concentrations  the  chemical  details  of  the  hydrogel  play  a 
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secondary  role  in  the  behavior.  Thus  autonomic  response 
depends  on  extrinsic  factors  including  the  reactant  concentra¬ 
tions,  the  structure  of  the  active-inactive  node,  and  tempera¬ 
ture.  For  example,  the  node  depth  is  limited  to  a  small  fraction 
of  the  overall  hydrogel  thickness  due  to  a  combination  of  factors 
including  overall  substrate  thickness  required  for  handling  and 
a  trade-off  between  ink  penetration  and  precision  feature 
printing.  This  constrains  the  swell-deswell  behavior  of  the 
printed  nodes,  and  thus  mechanical  swelling  behavior  was  not 
observed  in  our  composite  hydrogels.  Our  findings  suggest  that 
to  achieve  swelling  sufficient  for  two  active  nodes  to  commu¬ 
nicate  through  strain  distortions  of  an  intervening  inactive 
matrix,  the  minimum  penetration  depth  is  likely  to  be  ~40%  of 
the  total  thickness  of  the  sample.  Optimally,  one  would  create  a 
cylindrical  pillar  that  extends  the  entire  thickness  of  the 
sample.  The  reactive  printing  techniques  described  here  though 
would  not  be  suitable  for  such  a  design. 

We  demonstrate  that  temperature  exponentially  impacts  the 
oscillation  period  within  an  active  node.  Temperature  gradients 
direct  wave  propagation  through  stripes  or  arrays  of  active 
nodes.  The  shorter  waves  that  originate  at  a  higher  temperature 
eventually  dominate  over  those  from  lower  temperature  in  a 
manner  analogous  to  the  pacemaker  effect  seen  in  high 
frequency  nodes  in  close  proximity.'^  This  implies  that  inho¬ 
mogeneities  in  temperature  can  be  combined  with  mechanical 
swelling  in  a  way  similar  to  that  discussed  with  photo  illumi- 
nation^'’^^  to  create  directed  motion  and  transport  along  BZ 
reactive  strips. 

The  reactive  Ru  ink  described  here  can  be  used  to  build 
complex  patterns  of  active-inactive  composite  BZ  hydrogels. 
This  precise  control  of  spatial  patterning  enables  future  studies 
into  the  interaction  of  Ru  placement  with  other  known  triggers 
of  the  BZ  system,  such  as  light,  temperature,  and  mechanical 
strain."’'^'"  The  primary  limiting  factor  of  this  study  was  the 
trade-off  between  ink  penetration  and  pattern  precision.  Future 
alternative  Ru  reactant  and  polymer  chemistries  may  offer 
improved  binding  efficiency  and  solubility,  which  would  help 
mitigate  the  penetration  limitation.  The  merging  of  autonomic 
hydrogel  composites  with  the  rapid  maturation  of  additive 
manufacturing  tools,  including  inkjet  printers,  has  vast  poten¬ 
tial  due  to  the  ease  of  programing  and  creating  arbitrary  2D  and 
3D  structures.  In  the  future,  by  combining  reactive  printing  with 
different  solvents  and  substrates  that  enhance  mechanical 
swelling,  a  new  class  of  fast  fabricated,  autonomously  shape 
changing  devices  can  be  envisioned. 
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